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Abstract 

This paper considers a dynamic coverage problem for sensor networks that are 
sufficiently dense but not localized. Only a small fraction of sensors may be in 
an awake state at any given time. The goal is to find a decentralized protocol 
for establishing dynamic, sweeping barriers of awake-state sensors. Following 
Baryshnikov-Coffman-Kwak fiTj, we use network cyclic cellular automata to gen- 
erate waves. This paper gives a rigorous analysis of network-based cyclic cellular 
automata in the context of a system of narrow hallways and shows that waves of 
awake-state nodes turn corners and automatically solve pusuit/evasion-rype prob- 
lems without centralized coordination. 

As a corollary of this work, we unearth some interesting topological inter- 
pretations of features previously observed in cyclic cellular automata (CCA). By 
considering CCA over networks and completing to simplicial complexes, we in- 
duce dynamics on the higher-dimensional complex. In this setting, waves are seen 
to be generated by topological defects with a nontrivial degree ( or winding num- 
ber). The simplicial complex has the topological type of the underlying map of the 
workspace (a subset of the plane), and the resulting waves can be classified coho- 
mologically. This allows one to "program " pulses in the sensor network according 
to cohomology class. We give a realization theorem for such pulse waves. 

1 Introduction 

A wireless sensor network (WSN) consists of a collection of sensors networked via 
wireless communications, with every sensor being a device collecting data of the en- 
vironment with respect to one or more features, and returning with a signal l22l [Tl. 
Sensors can read, inter alia, temperature, pressure, sound, target presence, range, and 
identification. Current-generation smart sensors, increasingly smaller in size, can per- 
form data processing and computation, albeit with very limited memory and compu- 
tation capability. However, constrained by the locality of sensors' sensing function, 
networks of sensors have many more applications in monitoring larger domains than 
a single sensor. They can communicate with each other in the sense of transmitting 
and receiving signals, which allows local information to be collected and agregated 
globally. 
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Figure 1: Greenberg-Hastings model in square with blocks: black blocks are regions 
where no sensors are located, waves behave the same as with no blcks. 

A very common application of wireless sensor networks is intrusion-detection: the 
network monitors an area, reporting the existence of intruders when they are detected 
by at least one sensor. Video surveillance provides one such example. There is consid- 
erable activity in this field, focusing on different features and goals, "optimizing" net- 
works in various senses. One such aspect concerns coverage problems, which consider 
whether a domain is always fully covered by the union of sensing regions of sensors, 
static or mobile. Current approaches include methods from graph theory, computa- 
tional geometry, and algebraic topology fl5l[r7lfT4l l5l. Other aspects focus on provid- 
ing a specific degree of coverage, while keeping the connectivity of the network lEOl . 
What concerns us most in the present work is the minimization of energy consump- 
tion, keeping in mind that sensors are almost always battery-driven. One of the most 
intuitive ways is constructing a sleep-wake protocol for the network, allowing sensors 
to alternate between higher and lower energy cost states ||2) . 

This paper is motivated by recent work suggesting the use of cyclic cellular au- 
tomata (CCA) for intrusion-detection sensor networks J2). This inventive paper applied 
the Greenberg-Hastings automata on a two dimensional plane to generate "waves" of 
on-state sensors for intruder detection. Generally speaking, this specific automata as- 
signs to each sensor the state space Z n (the cyclic group on n elements), and the sen- 
sors update their state by advancing one automatically, except in state 0, in which case 
update to state 1 is induced by contact with at least one neighbor in state 1. This 
Greenberg-Hastings automata on lattices has been frequently investigated in the litera- 
ture as follows. Some rigorous statistical results has been proved for GHM with state 
space Z3 [7|. For general state space Z„, experiments has been carried out in 10, and 
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specific features or patterns that would keep and ergodic behaviors has been studied in 
(9] [10] |6) . A few authors have considered what happens on a graph as opposed to a 
lattice [16 1; this is the starting point for the application in fl2) to sensor networks. The 
main features and results of Baryshnikov-Coffman-Kwak include: 

1. The CCA runs on a random graph instead of a lattice, where nodes represent 
sensors and edges communication links. 

2. The network is completely non-localized and coordinate-free 

3. The CCA with random initial conditions generates the familiar spiral-like wave- 
fronts that sweep the whole domain with on-state sensors, giving a decentralized 
scheme for low-power dynamic barrier coverage. 

4. Parameters such as wavelengths are controllable. 

5. For planar domains with small obstacles, the wavefronts behave as if there are 
no obstacles at all (see Figure[T): waves propagate through, making the problem 
of undersampling ignorable. 

The present paper begins where J2) ends, by investigating what happens when this 
protocol is adapted to an indoor network where the geometry and topology is not that 
of an open plane (perhaps dotted with obstacles), but rather a system of fairly narrow 
hallways connected with a non-trivial topology: a "fat" planar graph. The contributions 
of this paper include the following: 

1 . We observe and then prove that wavefronts propagate through the hallways, turn- 
ing corners and branching off to side-corridors. 

2. We lift the CCA dynamics from the network to the higher-dimensional simplicial 
complex the network bounds. 

3. We detail a pursuit-evasion game within the domain and give sufficient condi- 
tions (in terms of the topological features of the system) for the pursuer to win. 

4. We show how wavefronts have well-defined cohomology classes and prove a 
realization theorem for which cohomology classes can be attained by the system. 

5. We identify what we believe is a novel type of global defect in CCA, generated 
by the topology of the domain as opposed to a local singularity. We show how 
such pulse solutions behave like solitons in the system. 

Most existing work on conserving energy for WSN focuses on distributed sleep- 
wake scheduling. For example, PEAS[21 1 provides a protocol by forcing a node who 
has an active neighbor to sleep for period according to exponential distribution. It is ro- 
bust against node failure, however, could not guarantee, or measure the coverage with 
the rapid change of active sensors. The CDSWS[ 18 1 protocol uses a clustering tech- 
nique to divide the sensors into multiple clusters, and selects a few sensors from each 
cluster to work, while maintaining nearly full coverage. ASCENT[3| allows sensors 
to measure their connectivity in the network in order to activate their neighbors based 
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on those measurements. But it never allows working sensors to go back to sleep again, 
which ends up consuming more energy as time goes by. Compared to those works, our 
protocol provides the "user" a chance to determine how much energy they would allow 
to be consumed, as balanced against how long it takes the system to detect evaders in 
the environment. The more energy it consumes, the less the expected time would be. 
Another advantage over the other protocols is it guarantees the failure of any evader 
following continuous path in the domain. Although our scheme requires synchroniza- 
tion ahead of time, and has not taken into account node and link failure yet, it provides 
a new approach to designing distributed sleep-wake WSN with energy constraints. 

The outline of our paper is as follows, ^provides a network protocol, along with 
simulations and observations. In fj3] we first introduce the topological tools used later, 
then classify the asymptotic behavior of the system, and detail a necessary and suffi- 
cient condition for the system to not converge to an all-wake state. Degree as a time 
invariant is first introduced. We also formally define the "evasion game" at the end of 
this section. <|4] verifies the system on the one-dimensional limiting space defined by 
the hallways. For proofs of normal hallways case as in the simulations, please refer to 
f|5] where we turn to topological tools. f|6]is a supplementary section, answering ques- 
tions such as what the dynamics would be if there is no "local defect", by paring degree 
and cohomology. Next comes Sf7] which briefly analyzes the system's feasibility even 
under link failures. Conclusion and comments are in <|8] 

2 Greenberg-Hastings Model and Simulation 
2.1 Cellular Automata 

A cellular automata, (CA), is a lattice-space and discrete-time dynamical system. 
Spatial coordinates are called nodes, and the dynamics generally take values in a finite 
alphabet A, with A = Z2 = {0, 1} being the most common choice. The dynamics 
are local, in that the update rule for a node is a function of its state and the states of 
its spatial neighbors. For example, in 1? lattice, the Von Neumann neighborhood of a 
node with coordinate is defined as the set of nodes attached to it, including itself, 
i.e., {(i ~ - 1), (i + + 1), An initial state (time t = 0) 

is selected by assigning a state for each node, typically at random. A new generation 
is created (advancing t by 1), according to some fixed rule universally that determines 
the new state of each node in terms of the current states of the node and its neighbor- 
hood. Typically, the updating rule is the same for each node and does not change over 
time, and is applied to the whole space simultaneously (but see asynchronous cellular 
automata [2| for one exception). 

This paper focuses exclusively on cyclic cellular automata(CCA). The alphabet is 
defined to be A — Z n = {0, 1, . . . , n — 1} under modular arithmetic. One denotes the 
(discrete) collection of nodes as X and denotes a state at time t as u t : X — > Z„. The 
updating scheme for a general CCA is increments states in Z„, assuming that some 
excitation threshold is exceeded. More specifically, u t+ i(x) = u t (x) + 1 if certain 
criteria concerning the states of the immediate neighbors of x, Af(x), are met. Such 
systems tend to cause periodic or cyclic behavior, spatially distributed and organizing 
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into waves. 

The Greenberg-Hastings Model (GHM) is a CCA first invented to study the spa- 
tial patterns in excitable media MTZj . In this model, special significance is given to a 
single state (state 0), interpreted as an excitation state. The update rule for GHM is as 
follows: 



This updating scheme is interpreted as the result of two mechanisms combined, 
excitation mechanism diffusion. It is therefore no surprise that there is a strong resem- 
blance between the behavior of GHM on the plane and solutions to reaction-diffusion 
PDEs on planar domains [11|, with both generating spiral-type waves. Given a fixed 
network, the states of the nodes will be uniquely determined by the initial state, because 
the GHM is a deterministic model. Denote by Q the evolution operator Q : — > Z„ . 

2.2 Observations 

Figure |2 illustrates the dynamics of the GHM on a specific indoor network. The net- 
work is built on narrow hallways modeled as a metric space with Euclidean metric; 
the neighborhood of a node is defined as nodes within distance r. We parameterize 
the system of 16250 nodes inside a 200 x 200 square with n = 20 and r — 1.5. The 
colors are representing states, with dark blue representing state 0. At time=0, it is in 
an unordered initial state. During the first several time steps, generally until time 25, 
the ratio of nodes with states grows, as a result of the fact that states grow steadily 
until they reach and wait for a stimulation from its neighborhood. At around time 45, 
spiral patterns become clearer visually, from top left, bottom and middle right. Those 
spiral "seeds" propagate waves along hallways. Wavefronts, consisting of the nodes 
with state 0, sweep through the domain, traveling "intelligently," turning corners, etc. 
When wavefronts coming from different directions meet, they annihilate. And after 
enough steps (about 250), wavefronts have finally covered the whole space. 

This protocol has some properties that make it ideal as a intruder detection sensor 
network. 

1. The system has tunable energy efficiency. Set state to be the waking state, with 
all other states 2, 3, . . . , n — 1 denoting sleep mode. Those sleep-mode sensors 
are doing nothing but advancing their states by 1 for every time step. This could 
be done with very low energy consumption because they only have to follow 
clock clicks with no computing, sensing, or transmitting. Intrusion-detection 
is performed by the wake-state nodes. After a sufficient time lapse, only a 
fraction (about l/n) of the total sensors will be in wake mode at any given time: 
the larger n, the less energy consumed. 

2. The wave length is generally fixed no matter which source it is from, assuming 
the nodes are uniformly distributed: it appears to depend linearly on n. For 
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Figure 2: Greenberg-Hastings Model for Narrow hallways space at time 0, 20, 45, 90, 
150, 200, 250, 350. 
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bigger n, longer wave length is generated, but seeds are generated with smaller 
probability and longer generation time. This makes a trade-off between energy 
consumption and system success. 

3. If we are given some specific sensors, say, with a big enough sensing radius e, 
then we will see that nodes in wavefronts form barriers, cutting the hallway into 
disconnected pieces. We note that the wavefronts efficiently sweep the corridors. 
Any intruder between two barriers has to follow the direction wavefronts propa- 
gate in order not to be detected immediately, but still is not able to survive to the 
end and will be detected by a coming wavefront in the opposite direction. 

3 Topological tools and dynamical features 
3.1 Topological tools 

Our goal of a decentralized protocol for integrating local data into a comprehensive 
understanding of the global system points to algebraic topology as an appropriate and 
useful toolset. In order to construct a topological object based on the network that pre- 
serves local information, we consider the flag complex of the network graph. Recall, 
a simplicial complex is a union of simplices obtained by gluing them together along 
faces of same dimension lfl3l . The flag complex (also known as a clique complex) of 
a network is the maximal simplicial complex with the network graph as 1 -skeleton. 

Definition 1. Given an undirected network graph G, with vertex set X and edge set 
E, the flag complex Cf, of G is the abstract simplicial complex whose k-simplices 
correspond to unordered (k + \)-tuples of vertices in V which are pairwise connected 
by edges in E. 

If metric information about the network is given, the Vietoris-Rips complex can 
also be built. 

Definition 2. Given a set of points X in a metric space and a fixed parameter r > 0, 
the Vietoris-Rips complex of X, lZ r (X), is the abstract simplicial complex whose k- 
simplices correspond to unordered (k + l)-tuples of points in X which are pairwise 
within distance r of each other. 

Compared to a flag complex, the (Vietoris-)Rips complex requires metric informa- 
tion about the space. But if the network on a metric space is built with communication 
radius r, which means two nodes are neighbors if and only if they are with in distance 
r, then the flag complex and Rips complex built on this network are exactly the same. 

Definition 3. For an abstract simplicial complex C whose O-simplices are located in a 
d- dimensional Euclidean space K d , the shadow of C in T>, S(C), is the union of convex 
hulls of each simplex in C. 
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3.2 Dynamical features 

We reprove certain results from the CCA literature [ 8 , 7 1 in the more general setting of 
network (as opposed to lattice) systems. Our perspective is that a CCA is a discrete- 
time network-based dynamical system. From observation, the interesting dynamical 
features associated with the GHM are time-periodic. We therefore focus our efforts on 
understanding time-periodic states. 

Definition 4. An orbit of a node x £ X under GHM with an initial state uq is the 
time-sequence of states (ut(x))£2. . A node x is said to be ^-periodic if its orbit 
satisfies Ut+K {x) = Ut (x)for some K > and all t. A node x is said to be eventually 
periodic if its orbit satisfies Ut+K{x) — Ut(x) for some K > and all sufficiently 
large t. 

Definition 5. A state u on a subgraph X' C X is continuous if for every pair of 
neighbors x, y £ X', \u(x) — u(y)\ < 1 (where, recall, all addition is in Z„J. 

Definition 6. A node x is subordinate to a neighbor y at time t if their states at that 
time satisfy Ut(y) = u t (x) + 1 (where, recall, all addition is in Z n ). 

Lemma 1. Subordinate nodes will remain continuous for all future time. 

Proof. It suffices to assume a subgraph consisting of a single edge with nodes x and y. 
Assume that \u t (x) — Ut(y)\ < 1. Consider the set S = {z £ X\ut+i(z) — u t (z)} C 
u t ' 1 (0). Depending on membership in S, 

!0 x, y £ S or x, y S 
1 x g S and y £ S (1) 
-1 x £ S and y g S 

By continuity, \u t (x) — u t (y)\ < 1, so |u t+1 (a;) — u t+1 (y)\ will exceed 1 only if 
u t (x) - u t (y) = 1, ((£/(«)) (a?) - (G{u))(y)) - (u(x) - u(y)) = 1 or it (at) - u(y) = 
— l,((Q(u))(x) — (Q(u))(y)) — (u(x) — u(y)) = —1. The first case is equivalent to 
u(x) — 1, u(y) — and x £ S,y £ S, which is impossible because y has neighbor 
x in state 1, and will not stay in state for the next step, thus not in S; the second 
case is the symmetric case which by the same argument is not possible either. Then 
\(Q(u))(x) — (G(u))(y))\ < 1, which makes Q(u) also continuous. □ 

Corollary 1. If a node x is subordinate to a neighbor y that is n-periodic at time t, 
then node x is n-periodic for all future time. Any node subordinate to an eventually 
periodic node is eventually periodic. 

Proof. Suppose x reaches for the first time (after t) at time to- By the scheme of 
G, u to (y) = 1. Therefore, all we need to prove is for any non-negative integer k, 
ut +kn+i{x) = 1. We have already proved u to +i(x) = 1 because it has a neighbor 
y in state 1 at that moment. Suppose the statement holds for a particular fco, i.e., 
ut +k n+i(x) = 1, then u to+{ko+1)n (x) = 0. But by periodicity, u to+{ko+1)n (y) = 1, 
thus u to +(k +i)n+i(x) = 1, which makes the statement hold at fco + 1. By induction, 
the statement holds for all k. □ 
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Figure 3: Counter example: non-continuous state for all time. 



Corollary 2. Continuity is forward-invariant: continuous states remains continuous 
in time. 

Proof. According to Corollary [T] two neighbors that are subordination will remain 
continuous. For one step forward, two neighbors that are of the same state will either 
remain the same state, or be offset by state 1, which means subordination, thus also 
continuous. □ 

However, it is not necessarily the case that all initial conditions converge to a con- 
tinuous state (even in a connected compact network). See, for example, FigureO every 
node has period n, and the two nodes on the right end have states always differed by 3. 
Thus this is never a continuous network. 

The following definition is a network-theoretic version of the lattice-based ana- 
logue from, e.g., (8). 

Definition 7. We call a formal linear combination a of edges on = [a,i,bi],i = 
1, . . .,K a cycle if the boundary of a, da = (bj — ai) is 0. A cycle is called 

a loop if hi — ai + 1/or i = 1, . . . , K — 1, and bx = &1- 

As a remark, a loop is a cycle, and a cycle is the sum of one or more loops. We 
also remark that the set of cycles Z has the structure of an abelian group: one can add 
cycles and scale them by (integer) coefficients. 

Definition 8. A state u : X — > Z„ has a seed if there is a loop Xi^o^^'^+i] 
(xq — Xk), for which u{xi) = i mod n. 

By definition, the length if of a loop that makes a seed has to be a nonzero multiple 
of n, because u(xo) — u(xk), K = mod n. Every node on a seed has period n, 
because it will have a neighbor of state 1 on the seed when it reaches state 0. 

Lemma 2. If an initial condition uq on a connected compact network X contains at 
least one seed, then all nodes are eventually periodic. 

Proof. Let the set of n-periodic nodes in X be P t , and state at time t be ut- Suppose the 
loop X^q 1 [ Xi > } x o = x k makes a seed in initial condition, then Pq is nonempty, 
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Figure 4: A seed (left) and a defect (right) for n = 8 on cycles in light red. 



with (xi)g _1 as a subset. Pt is non-descending with respect to time t, Pq C Pi C ■ ■ • C 
P t C Pt+i C . . . . For a node x that is not n-periodic that has at least one neighbor 
that is n-periodic at time to, if x never gets to be n-periodic, it means for any t positive, 
there exists some s > t such that u s+ i(x) ^ u s (x) + 1 mod n. It will induce that 
u s+ i(x) = u s (x) = 0, which is saying x gets to stay in state for a while from time to 
time. But the neighbors of x that are periodic n are advancing their states by 1 at every 
time step, this will make the face difference between them and x bigger and bigger until 
it reaches 1 mod n. When the such a offset by 1 appears, a subordination between x 
and its periodic n neighbor is built up, which makes x periodic ever since as a result of 
corollaryQ] Therefore every node which as at least one neighbor that is periodic n will 
be periodic n after a finite amount of time (no longer than n). By the above argument 
and the fact that the network is connected, for any t, Pt C P t+n . On the other hand 
side, since X is compact, there exist a time T, such that 



We see in the above arguments that a loop that makes a seed at one moment will 
support a seed forever with the dynamics. The key feature that is invariant under the 
dynamics is the concept of "winding number", which records how many rounds it goes 
through while chasing continuously on a loop. We will define this as degree and extend 
the concept to all cycles. 

Definition 9. For a given network X and a state u £ Z„ , ifu is continuous on a cycle 
a = J2i=i[ a i> bi], then the degree of u on this cycle is defined as 



|J Pt = Pt 



Therefore, the whole system is in ?i-periodic state since time T. 



□ 



fc-i 




i=0 
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where the summands are forced to be — 1, 0, or 1, and the sum is ordinary addition ( not 
mod n). 

Definition 10. We call a cycle a — in the network X a defect for some 

state u € Z„ if the degree ofu on this cycle is nonzero. 

An example of a defect is as in Figure[4] The concept of a defect is a generalization 
of a seed, in the sense that it has nonzero degree. The term "degree" defined here is 
consistant with the use of degree in topology, which is a homotopy invariant [ 13 1. Here, 
it is the discrete version of "winding number" for continuous self-maps of the circle S 1 
[7 1, describing how many times it wraps around with direction. Similar to Lemma 
5 in 0, we will prove the R 2 version instead of the lattice 1? version, presenting a 
necessary and sufficient condition for a continuous system not dying out. 

Lemma 3. For two cycles a and f3, if a state u is continuous on both cycles, then it is 
also continuous on their sum a + f3, and deg(u, a + /?) = deg(u, a) + deg{u, j3). 

Proof. Let a = J2i=i b i\ and I 3 = Yn=i \ c h <k], then 

deg(u, a + j3) (2) 
fc-i fc-i 
=l/n(J2(u(h) - u{a,)) + ^O(di) - u(ci))) (3) 

i=0 i=0 
fc-1 fc-1 



=l/n^j(u(&i) - u( ai )) + l/n^2(u(di) - u(c t )) (4) 

t=0 i=0 

=deg(u, a) + deg(u, f3) (5) 

□ 

Lemma 4. For a cycle a and a continuous state u, the degree of u on this cycle is 
invariant under the GHM updating rule Q, i.e., 

fc-i fc-i 

- (S(u))( Xi )) = l/n^2(u(x i+1 ) - u(xi)) (6) 

i=0 i=0 

Proof. We first prove that degree on a loop X^o 1 i Xi > x «+i] ; ^o = %k is invariant. As 
before, equation Q] holds for every pair of neighbors Xi+\ and Xi. Since the number 
of pairs (xi+i,Xi) with Xi E S, Xi+i ^ S is the same as the number of pairs with 
Xi £ S,x i+ i E 5,the summation of ((g(u))(x i+1 )-(g(u))(xi))-(u(x i+1 )-u(xi)) 
is 0, which makes Equation|6]hold. Since every cycle is the sum of one or more loops, 
and degree is additive by Lemma[3] then it is also invariant on cycles. □ 

For a state u on a loop that forms a defect, if the loop bounds a region V in R 2 
that belongs to T>, we can discuss the continuity of the subnetwork in V . If the subnet- 
work in V is sufficiently dense (e.g., 2-complex has shadow containing V), we observe 
that the subnetwork could never reach continuity, with at least one singularity (a dis- 
continuity) forced, as in Figure This could be understood intuitively as a discrete 
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Figure 5: In the region bounded by a defect, the state is discontinuous. 



version of the theorem in complex analysis, which says a holomorphic function on a 
domain always has integration on the boundary. It would also contradict the fact that 
a continuous map from a contractible space to S 1 has degree restricted on any loop. 

Lemma 5. Consider a state u on X with n > 3, and a loop I = $^_n 1 [zEi, ^i+iL — 
xk in X. If the loop I is null homologous in the 2-complex built on X, and u on I makes 
a defect, then u is discontinuous on X. For n < 3 any state on X is continuous. 

Proof. Since I is null homologous in 2-complex built on X, I = d(3 where (3 = 
J2iLi A i s a 2-chain in the 2-complex and are 2-simplices. I can be deformed 
to a single 2-simplex through a sequence of homologous loops in X, Y^h=\ dfiiiJ = 
1, . . . , m — 1, while the successive two loops only differ by the boundary of one 2- 
simplex. Suppose u is continuous on X, such operation could not change the degree at 
all, since at most three of the summands u(xi + i) — u[xi) has been changed value up to 
1 . Thus the summation is at most changed by 3, which makes the degree changed by at 
most 3/n, which has to be invariant when n > 3. Thus the nonzero degree remains the 
same for the sequence of homologous loops, which can not be true because the degree 
on the boundary of a single 2-simplex has to be zero. Therefore the state on X could 
not be continuous. 

It is trivial to see the continuity when n < 3, because any two elements in Z„ differ 
by at most 1 . □ 

Theorem 1. For a continuous state u on a connected compact network X, the system 
eventually turns to all-0 state (die out) if and only if u does not contain a defect. 
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Proof. Suppose u contains a defect on cycle a = J2i=i [ a ii By l emma IU the degree 
is invariant under Q, so it will never be 0, thus the system will never turn to all-0 state. 

For the converse, we need to show for a continuous state u not dying out eventually, 
it has to contain a defect at the beginning. Firstly, it is obvious that after long enough 
time, in such system w t , every node in state i must have a neighbor in state i + 1, for 
all i 0, since Ut = G z (ut~i). So if we start from a node xq, such that u t (xo) = 1, 
we can find a neighbor of xq, x%, such that u(xi) = 2. Following the process, we 
get a sequence of nodes Xo,X\, . . . ,x n -i, such that Xi and Xj+i are neighbors and 
Ut{xi) = i + 1 mod n. If from every state node, a state 1 node could be reached 
by jumping along neighbors which are in state 0, then following the process, we will 
finally reach a node has been visited before, as X is compact. In this way, we have 
obtained a defect in Ut, which is also a defect in u, by lemma [4] and the fact that u is 
continuous. 

To see that a node with state 1 could always be reached from a node with state 
by jumping along neighbors in state 0, all we need to prove is there is no such set A of 
nodes with state 0, that their neighbors not in A could only be in state n — 1. If such 
A exists in u t , then there is a proper subset of A with state 0, and state n — 1 on the 
complement in Ut-%. Following these procedure, we should finally obtain a set Aq of 
state nodes, each has at least one neighbor with state n—\ and other neighbors with 
state in u s . Then in u s _i, nodes in Ao have to be in state n — 1 (by continuity), and 
their neighbors not in Ao must all be in state n— 2, and for n steps back, in u s - n , nodes 
in Ag have to be in state 0, and their neighbors not in Ao must all be in state n — 1. But 
such a u s _„ could not produce u s _ Jl+ i under Q, because those state nodes have no 
neighbors in state 1. Therefore such a set A does not exist, which makes the statement 
in the beginning true. □ 

Since degree is invariant under the update rule Q, Theorem[T]can be interpreted as 
saying that a continuous state dies out eventually if and only if it is cohomologically 
trivial (see f|6]for details on how to define the cohomology class of a state). 

Theorem 2. There exists a directed subgraph J- of the network that is a spanning 
forest rooted at seeds, with directed edges in the direction of subordination. 

Proof. Every node that is not originally a seed node will become periodic by building 
up a subordination with some periodic neighbor. For every non-seed node, choose 
one from its neighbors via subordination and use a directed edge with itself as head 
to represent the relationship. This forms a directed subgraph of the network. From 
any non-seed node, following those directed edges with inverse directions, it has to 
end in a seed node, because it is a compact network. We argue that the subgraph is a 
tree because it contains no loop: if it did contain a loop, then the loop is comprised of 
non-seed nodes, but for any directed edge, the head node becomes periodic later than 
the tail node, which is a contradiction with being in a loop. And furthermore, we can 
treat the forests as rooted at seeds, which makes every edge in the direction that goes 
deeper in a branch to the leaves. Such structure gives the nodes a hierarchy, and since 
for every edge, the two end nodes have states offset by 1, we can induce the state after 
the system reaches equilibrium. □ 
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Definition 11. The depth of a node in a tree is the number of hops between the node 
an the root of the tree. 

According to the above proof, we have made a point in that the growth of the forest 
is at most one level at a time, which means in every time step, there could not be more 
than one node from a same branch that becomes subordinated. 

Starting from a uniformly randomly generated initial condition (a reasonable if ide- 
alized statistical model) the system is not guaranteed to converge to a periodic system, 
not including all-zero states. One sufficient condition is the existence of a seed, which 
we prove to be of high probability with certain reasonable assumptions (Lemma [6}. 
It is possible that the system became messy with no wavefront observable (too many 
seeds all around in the space, for instance). We would require those nodes that are far 
away (in the hop-metric) from the defects to be in state at one moment (in our case, 
larger than the number of states is already enough). This assumption is proved later to 
be of high probability (Lemma|7]l. Under the above two assumptions, continuity in the 
acquired region will be guaranteed. Therefore, from now on, we limit discussion to the 
region far away from defects. 

Lemma 6. For a fixed uniformly sampled network with communication radius r and 
fixed n on a domain consisting affixed narrow hallways, the probability of at least 
one seed existing in the initial condition generated according to uniform distribution 
approaches 1 as the number of nodes grows. 

Proof. Divide the space into square shaped pieces T>i indexed by /, with side length 
smaller or equal to r/^/2. As the network size \X\ approaches infinity, the probability 
of there to be no less than n nodes in each T>i approaches 1 . For a T>i with n or more 
nodes, the sub-network in this subdomain makes a complete graph. Therefore, there 
is no seed in the initial condition, if and only if the nodes do not cover all the states, 
which means there is at least one state missing in the initial condition. Thus 

P(no seed in initial condition in T>i with m, nodes) 

^ n(n- l) m * 

= n(l - l/n) mi 

(7) 

and 

P(no seed in initial condition in T>) 

< J^J P(no seed in initial condition in T>i with m, nodes) 

i=l 
\I\ 

< Y[n(l - l/n) m * 

i=l 

= nl J l(l - 

(8) 
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which approaches as \X\ approaches infinity. 



□ 



Lemma 7. Starting with a fixed network and uniformly distributed initial conditions, 
with probability approaching 1 as the state number n grows, nodes with hop distance 
to all defects bigger than 2n will turn to state after 2n — 2 time steps. 

Proof. Suppose there is no state 1 node in u n -i in the region n hops away from any 
defect, then there could be no state 1 or 2 node in u n in the region n + 1 hops away 
from any defect, and with the same reason, there could only be state node in u 2n -2 
in the region 2n hops away from any defect. Therefore the probability that every node 
at least 2n hops away from defects are state in u-m-i is no smaller than the that of 
no state 1 node at least n hops away from defects in u n -\. 

Now suppose there is a node x at least n hops away from any defect, and u„_i (xi ) = 
1. Such x\ must have at least one neighbor of state 2, named x 2 , otherwise in one step 
before, it would not be able to update from to 1. Via the same argument, there exists 
a sequence of nodes: Xj, j = 1, 2, . . . , n — 1, such that Xj and Xj+i are neighbors, and 
u n -i(xj) = j. For one step ago, u n -2(xj) = j — 1 for j ^ n, and two steps ago, 
u n s(xj) = j — 2 for j = 2, . . . , n — 1 and u n -z{x\) G {0, n — 1}. Following such 
argument, back at time 0, Uo(xj) G {0, n — 1, . . . ,j + 1} for j = 1, . . . , n — 2 and 
uo(x n -i) = with at least one neighbor of state 1. 

Let Ij = {0, n — 1, . . . , j + 1}. For a fixed node x, 

P(at least one of x's neighbor have a state in Ij at time 0) 

(9) 

where is the number of neighbors of node x. Suppose N is a universal upper 

bound on |A/"(a;)|, then 

P (u n -i {x) = 1 for some x at least n hops away from any seed) 

n-1 

<\X\H(l-(l-j/n)^) 

3=1 
n-1 

<\X\H(l-(l-j/nf) 

3=1 

< \X\(l - (l/2)^)("- 1 )/ 2 

(10) 

which approaches as n approaches infinity. □ 

For example, in a 40000 nodes network, where every node could have up to 6 
neighbors and n = 20, the probability of a seed is bounded below by 0.9656, which 
validates the observation in previous simulation. As per the above two lemmas, we will 
always assume that at least one seed exist in initial condition, and the nodes at least 2n 
hops away from any defect will turn to state after 2n — 2 step. These two assumptions 
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guarantee not only the system not dying out (turn into an all-zero state), but also the 
continuity of the system in acquired region, with the following corollary. 

3.3 Evasion Game 

We propose a sensor-network based "Evasion Game" formally, and then use the model 
to verify the system: why the wavefronts sweep the entire domain; how to interpret the 
phenomenon that wavefronts are dividing their neighborhoods and how to prove that 
an intruder will always fail to evading detection; what are the parameters that control 
the system and how they are changing the behaviors of those wavefronts. 

Definition 12. Let the domain where the evader and sensors are located be denoted 
V C R 2 , and the sensor network X. For each sensor x € X, its coverage is a subset 
U x C V. Denote by X(t) the set of sensors in wake-state (0) at time t. We define 
the Evasion Game as follows: the strategy for the pursuer is to control the network 
following GHM, and the strategy for the evader is to pick a moment to to come into the 
domain, and follow a continuous path in V: f : [to, oo) — > V. The pursuer wins iff 

3 t € [to, oo), such that 

fir) € (J U x . 

xeX([ri) 

Otherwise, the evader wins. 

We note that the only requirement on the evader is its trajectory be continuous: 
there are no constraints on the velocity or acceleration. Even with such minimal con- 
straints, the evader is not able to win. 

4 Limiting case with 1-d hallways 

We begin our analysis with the limiting case when every hallway is sufficiently narrow 
compared to the walls, so that the domain V can be approximated as a (topologically 
equivalent) one-dimensional space. We assume those sensors are located in T> with 
each node having a coverage which is a one dimensional convex set around itself, and 
the convex hall of two neighbors is covered by the union of their coverage regions. We 
also assume the union of convex hulls of neighbors (subspace of T>) is good enough 
to cover V, in which case the whole space is covered when every sensor is turned on. 
If we run GHM on this network, with at least one seed in initial condition, then every 
evader (not near the seeds) loses the evasion game. 

Theorem 3. For GHM on network X with communication distance r in a compact 
and connected 1-d complex V, if the initial condition contains at least one seed, and 
there exists a subnetwork X' covering a sub-domain V , such that the state on X' is 
eventually continuous and contains no defect, an evader will always lose the evasion 
game on V. 
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Proof. For any time to when the evader comes into the domain, consider the product 
space V x [to, oo) with the second coordinate representing time. Treat the coverage 
of the sensors also as a subspace P c of V x [to, oo), which is 

oo 

|J |J U x x[t,t + 1) . 

t=fto],t£Zx£X'(t) 

Let p be the projection map: p : V x [to, oo) — > T>',p(a, t) = a. Then p : P c — > V 
is onto, because V is fully covered when every node in X' on. If we could prove that 
there exists a subspace in P c homeomorphic to V, with map p as homeomorphism, 
then V x [to, oo)\P c contains no continuous path from top V x {to} to bottom V x 
{T} for T big enough, because they are dual to each other. Therefore no evader could 
survive forever. The construction of the subspace in P c is as follows in Lemma [8j 
below. □ 

As a remark, a good example for the state on subnetwork X' is eventually contin- 
uous and contains no defect is to let it be all-0 state at a moment, which is observed 
most of the time in simulations. 

Lemma 8. Under the conditions of Theorem\3\ there exists a subspace S £ P c , such 
that p induces a homeomorphism from S to T>' . 

Proof. First, reduce to a subnetwork X" of X' such that within X" the convex hulls 
of neighbors is still enough to cover V, but any two distinct convex hulls intersect in 
at most one node. We then construct S inductively from the empty set as follows: 

1. Select an integer time t which is no earlier than to big enough, such that every 
node in X" has already been periodic for a long enough time. Pick a node 
x £ X"{t) and add (x, t) to S. 

2. For any neighbor of x in X" , say y, there exists a continuous path lying in P c , 
between (x, t) and (y, t y ), where y £ X"(t y ) and |t — t y \ < 1 (t y is an integer 
time), which is mapped homeomorphically to the convex hull of x and y in V , 
because continuity holds on edge [x, y], and U x X [t, t + 1) U U y x [t y , t y + 1) 
is a path connected set. For any x's neighbors y that has not been visited, add 
(y, t y ) with the continuous paths between (x, t) and (y, t y ) to S. 

3. Repeat step 2 for every newly visited node, until every node in a connected com- 
ponent of X" has been visited. 

Such procedure could not be realized only if there is a cycle in X" , such that the 
continuous lift of the path to P c is not a loop, which means the state restricted on the 
loop is a defect. However there is no defect in ut (X") when t is big enough, by Lemma 
2] Therefore the procedure is well-defined. Start above procedures until every node in 
X" has been visited. 

Such an S is mapped onto V by p, because every convex hall of two neighbors, say 
x and y, is mapped onto from the path between (x, t x ) and (y, t y ). The restriction of p 
to S is also injective, because every node and edge is only visited once, and p restricted 
on every continuous path between (x, t x ) and (y, t y ) is homeomorphism. 
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Figure 6: Product space and S: space V on top, time increases from top to bottom. 
State space is Z4, with white for state 0, light green for state 3, red for state 2, and dark 
blue for state 1 . Gray curve represent S. S = V by p. Continuous curves connecting 
top and bottom without intersecting S do not exist. 

The only thing left to be proved is that there exist a continuous inverse of p\g. Let / 
be a map from V to S, such that / maps every node x in X" to (x, t x ) in S, and maps 
every edge between node x and y to the continuous path between (x, t x ) and (y, t y ). 
Such / is an inverse of p\g, and is continuous: for a point in V that is not a node, it's 
covered by a convex hall of two neighboring nodes in X" , thus its small neighborhood 
maps to the lift of the the convex hall in S continuously; for a node point x in X", 
its neighborhood maps to a neighborhood of the lift (x, t x ) homeomorphically, by the 
procedure of constructing S. Thus / is an continuous inverse of p restricted on S, thus 
p induces a homeomorphism from S to V . □ 

5 Main theorem and proofs 
5.1 Assumptions 

The main theorem for this paper, Theorem|7] shows that any evader in the evasion game 
on a narrow hallway space Pel 2 will lose, given the appropriate assumptions about 
the density of the network X and the initial condition. Specifically, we assume: 

1. the projection from Rips complex lZ r (X) to T> preserves homotopy type; 

2. each sensor x £ X covers a convex set U x C T> around its location; 

3. the convex hull of sensors that are pairwise neighbors is covered by the union of 
coverage of those sensors; 
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4. there is at least one seed in the initial condition. 

According to a theorem of [4 1 , we will be able to build the correspondence between 
the Rips complex TZ r (X) of a planar point set and its shadow S(7Z r (X)) in K 2 . 

Theorem 4. [ IS!}] For any set of points in R 2 , ni(lZ r (X)) -t m(S(U r (X))) is an 
isomorphism. 

Definition 13. A local hole in the Rips complex is a non zero element of m(TZ r (X)) 
that has trivial projection in iri(T>). 

In the sense of local holes, Theorem|4]is saying that in our case, the Rips complex 
lZ r (X) has no local hole if and only if its shadow S(lZ r (X)) has no local hole. 

Another useful fact is that with very high probability, when the network is dense 
enough, the Rips complex 1Z r (X) has no local holes fl9l - Therefore, if with enough 
sensors uniformly distributed in the domain and with high probability, the Rips com- 
plex lZ r (X), and its shadow S(lZ r (X)) both have no local hole. 

5.2 Wave propogation 

Definition 14. A boundary path along a boundary component, is defined as a simple 
path such that every node on the path has a coverage that intersects with the corre- 
sponding boundary, and the intersection of the coverage of every two neighbors, x and 
y on the path, also intersects the boundary nontrivially. A boundary of a network X, 
dX on T> is a collection of boundary paths, one with each component of dT>. Refer to 
Figure\Sybr illustration. 

Definition 15. For positive integer set A, define the depth A nodes, Xa, as the set of 

all the nodes that are with depth k € A in the directed forest J- built on the network. 

Definition 16. A connected sub network X' of X makes a barrier, if there exist a piece 
of hallway T>, which intersects dT> at &D, and the composition d o i» : Hi(T>, dT>) — > 
H (dV) ofu ■ Hi(V,&D) ->■ Hi(V,&D) and d : H\(V,&D) -> H (dV) is an 
injection, such that X"s coverage contains at least one element in a nonzero class of 
H\(T>, &D). In other words, it covers a region that divides the hallway locally and 
transversally as in Figure^ 

Theorem 5. Let X be a connected and compact network on a narrow hallway space 
T>, running under GHM, whose initial condition contains at least one seed; if there is 
a moment that the subnetwork X' in a sub domain V is continuous and contains no 
defect, and boundary paths dX' exists, then if at time t, there is a wavefront that makes 
a barrier, which is not supported on any end leaves of the forrest J-, then there would 
be a wavefront also makes a barrier at time t + 1. 

Proof. If there is a wavefront of nodes with depth k at that makes a barrier, then we 
want to prove that a wavefront of nodes with depth k + 1 exist which also makes a 
barrier. Let A denote the subcomplex on subnetwork X<k+t> an d let B denote the 
subcomplex on subnetwork X>k+i- Then An B is precisely the subcomplex with 
nodes with depth k + 1. On the other hand, A U B is the whole complex, because for 
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Figure 7: A piece of hallway with a subnetwork that is a barrier. 



every simplex in the whole complex, their vertices are pairwise neighbors, so by con- 
tinuity of states on X' (by the assumption that there is a moment that the subnetwork 
X' is in state zero), their depths could only differ at most by 1, by Corollary [2] which 
means the simplex is either in A or in B. The Mayer- Vietoris sequence for A and B 
gives: 



Let [a] € Hi (A, dX'), \fi] <E H x (B, dX'), where a and /3 are both connecting bound- 
ary nodes of different sides. Such an a exists because of the existence of a previous 
wavefront of depth k, and j3 exists because the network is sufficiently dense in V, and 
a is not supported on any end leaves of J 7 . Then V* ([<*]> [P]) = (if not, let j3 be of 
opposite orientation) in H\{A U B, dX'), because first homology of A U B is trivial. 
Therefore ^([a], [0]) and ^([0], [/?]) are homologous. Thus ([«],[/?]) <G kerip. By 
exactness, ker-0 = im <fi, thus there exist a 7, such that ^([7]) = ([a], [(3]). As <\> is in- 
duced by inclusion maps, 7 has to be a path connecting boundaries of two sides, which 
is a wavefront of depth k + 1, cf. Figure [8] Therefore, by induction, barrier-inducing 
wavefronts of every depth exist. □ 

The above wave propagation theorem presents how waves travel along the hall- 
ways, but did not mention the generation of waves. The following proposition would 
explain how a first wave is generated under same assumptions. 

Proposition 1. With the assumptions of Theorem^ at least one wavefront that is a 
barrier and intersects with boundary paths on both sides must be generated. 



H^AnB^X') 




H 1 (A,dX')®H 1 (B,dX') 




(11) 
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Figure 8: Wavefronts propagation: light green nodes are with depth k — 1, red nodes 
are with depth k — 1, dark blue nodes are with depth k + 1; Boundaries of the domain 
are covered by boundary paths. 

Proof. Let X< k be the set of nodes with depth less than or equal to k. Then there is a 
filtration of Rips complexes: 

K r {X {0} ) c ft r (X<i) c ■•• cK r (X< k ) c ...TZ r {X) 

Since they grow by attaching nodes within communication distance as k increases, 
therefore, as TZ r (X) is connected, there has to be a fco, such that 7?., (X<fc) are all 
connected for k > k^. For two boundary nodes of Xif- \, if they belong to boundary 
paths near different boundaries, since they are connected, and for the same argument 
from Theorem |5]by using the Mayer- Vietoris sequence, they are connected by a path 
with nodes from X^ g y This path generates a wavefront that is a barrier. □ 

The above results not only explain why the evader has to lose the evasion game, but 
also explain the behaviors of the wavefronts seen in simulations. After the first several 
steps, the nodes far away from the seed are all turned on, until wavefronts generated 
by the seeds reach them. The movements of wavefronts are verified to be moving 
away from seeds, and they provide locally separating barriers, as observed. Another 
significant property we observe from simulations is that the wavefronts make turns 
when reaching a corner, as shown in Figure [9] This reminds again that the behavior of 
the system does only depend on topology, not geometry, of the underlying space. 

5.3 Main theorem 

For now, we will start argue that under certain conditions, evader will always lose the 
evasion game. 



21 



Figure 9: corner of a hallway: with state space Z4, white for state 0, light green for state 
3, red for state 2, and dark blue for state 1. The outer side boundary path have more 
nodes than the inner boundary path, but more nodes stay in the same states: four in 
light green and three in red. Thus the wavefronts propagate from vertical to horizontal. 

Definition 17. For x G X, define the open star (JJ^ of x, U x C TZ r (X), as the union 
of x and all open simplices with x as one vertex. 

Lemma 9. Let a be a d-simplex in lZ r (X). If there is a continuous function f : 
[ti,t 2 ) ~~ * S(a), such that f(t) ^ Uiex(t)nir^' Vt ^ [ti,t 2 ), then there exists a 
continuous function f : [t±, t 2 ) — > a, such that f(t) (£ \J xe x(t)na Ux, Vt G [ti, t 2 ). 

Proof. For a 2-simplex a = [xo, x\,x 2 \. There exist a homotopy equivalence h from 
S(a) to itself, such that the interior of f] i=0 1 2 U Xi H S(a) is mapped onto the interior 

of S(a), which is f) i=0 1 2 U Xi , and the inverse image of every open edge S([xi,Xj)) 

belongs to U Xi n U x . n S(<r). Therefore, we can construct / as ho /, with the property 

that/i-^n^E^) C riieA U *<> which induces /(i) $ \J xe x(t)nv &x, Vt G [ti,t 2 ). 

□ 

Lemma 10. For the Rips complex lZ r (X), if there exists a continuous function f : 
[to, 00) — > S(lZ r (X)), such that f(t) ^ Usex(t) Ux, Vt G [to, 00), f/zen there exists 
a continuous function f : [to, 00) — > 1Z r (X), such that f(t) ^ V} X £X{t) Ux, Vt G 
[t , 00). 

Proof. The 2-complex C 2 (X) as a sub complex, has the same shadow as the Rips 
complex 1Z r (X). Lift the path / from the shadow S(C 2 (X)) to the complex C 2 (X), 
then apply lemma|9]on every simplex it goes through. This will give a lift / : [to, 00) — > 
C 2 (X) c ft r (X), such that f(t) i Uex(t) Vt G [t , 00). □ 
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Theorem 6. For a network X with coverage regions hi and Rips complex lZ r (X) with 
the shadow the whole T>. Then if the state on X is eventually continuous, and if there 
is no defect in the initial condition, then the evader loses the evasion game. 

Proof. Suppose there is a continuous path / for the evader to follow in order to win the 
evasion game, / : [to, oo) — > S(lZ r (X)) = T>, then by LemmafTOl there exists a lift of 
/, / : [to, oo) — > TZ r {X), such that following /, the evader could win the evasion game 
with coverage regions {U x \x € X}. Furthermore, since f(t) £ Usexft) &x> Vt G 
[0, oo), and by the fact that a 1 -simplex is covered by a subset of sensors that covers 
the simplex containing it, we can construct a continuous path /' that travels only on 
the 1-skeleton of lZ r (X) and still is safe, never being detected. However, by the same 
argument as in Theorem[3] since there is no defect in initial condition, such a strategy 
does not exist: any such evader would lose the game. □ 

If may not be the case that S(TZ r (X)) D T>. Our approach for solving this problem 
is by adding sensors to the network without changing the coverage, but enlarge the Rips 
complex such that it projects onto the whole domain. 

Lemma 11. If a boundary path exists within distance y3/2r to each boundary com- 
ponent of dT>, then there is a new sensor network X by adding sensors to X, with the 
same coverage at every moment, such that the shadow oflZ r (X) is T>. 

Proof. For every node x in the boundary path, add a node x' in U x n &D to the new 
network X, and for every edge on the path [x,y], add a node z' in U x fl U y n dT> 
to X' . For a quadrangle with vertices x,y,x',y', it is covered by union of U x and 
U y . Let x' and z' have same coverage and states as x, and y' has the same as y after 
equilibrium, then the coverage of X is exactly the same as that of X at every moment. 
Another property worth noticing is lZ r (X) now has its shadow same as T>, because 
'}, [y, z'], [z', y'}, [y\ y) are all 1-simplices in lZ r (X), which makes 
the shadow exactly T>. □ 

Theorem 7 (Main Theorem). With the existence of boundary paths within distance 
\/3/2r to boundary of hallways, the evader will always lose the evasion game in the 
sub domain T>' C T> on which the states is eventually continuous and contains no 
defects. 

Proof. By Theorem[6]and LemmafTTI □ 

6 Controlling the Cohomology 

We have observed in simulation that sometimes there is no "local defect" continuously 
generating wavefronts, but the system still reaches a nonzero equilibrium, with the 
remaining wavefronts propagating along hallways in periodic way. This phenomenon 
contributes to the existence of a "global defect", which differs from the "local defect" 
in that the cycle on which the defect is supported is in a non-zero class in first homology 
of the Rips complex lZ r (X), instead of a trivial one. Such an equilibrium presents a 
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Figure 10: add new nodes x',y', z' to the network, with x' and z' have same state and 
coverage as x, and y' has same state and coverage as y 

much higher portion of nodes in state than those with local defects. We will try to 
manually generate such patterns in GHM by turning off local defects. 

Such protocol is not energy efficient unless we shift state to sleep state as follows: 
the new interpolation lets state 1 to be waking state, state 2 to be broadcasting state, 
and state 3 till to be sleeping state. Then most of the nodes will be sleeping after they 
are eventually periodic. 

Recall from Definition^ the degree (or winding number) of a continuous state on 
a cycle is an index measuring how man times the states cycle through the alphabet on 
the cycle. Therefore, after local defects are turned off by breaking the links between 
state and state 1 nodes, the degree of a cycle which makes a nonzero class in first 
homology of the Rips complex lZ r (X) is determined by the number of wavefronts 
already generated and their directions of propagation. In other words, degree for all 
cycles is determined absolutely by local defects' location and the number of wavefronts 
they have sent out in the hallways. Note that the degree is invariant in time for a 
continuous state. Thus counting the degree for a cycle after defects are turned off is not 
a difficult problem: following the direction of this cycle, the number of wavefronts in 
the same direction minus the number of wavefronts in the opposite direction determines 
the degree. 

Definition 18. For a class [a] in Hi(1Z r (X)) (or Hi(C<z{X)), the first homology of the 
2-complex), define the degree of a continuous state u on [a] as the degree ofu on cycle 
a. If the projection ir : lZ r (X) — > T> induces an isomorphism : H\{TZ r {X)) — > 
Hi(T>), then define the degree ofu on ir*([a\) as the degree ofu on a. 

As a remark, the degree of a generator in first homology is well-defined, if for 
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homologous cycles a and /3, degree of u restricted on both are the same. By Lemma 
[5] degree of u restricted on a — (3, which is a null-homologous cycle, has to be zero. 
Thus, the degree on a and f3 have to be the same. By abuse of notation, deg(u, a) will 
be used for a as a first homology class in either lZ r (X) or V. 

Definition 19. Let Cont(X) represent the set of continuous states on X. Define a 
cohomologizing map h : Cont(X) — > H 1 (lZ r (X)) — Hom(Hi(lZ r (X)),1,), such 
that h(u)([a\) — deg(u, [a]). 

As a remark, the first cohomology H 1 (lZ r (X)) defined here is a simplicial coho- 
mology. It is torsion free and therefore can be treated as Horn(Hi(TZ r (X)), Z). 

Definition 20. A single wave is a continuous state u on X, such that (1) there exists a 
barrier on which u is supported, and (2) there exists a cycle a on which the degree of 
u is 1. 

By the definition of a wave, the degree is zero on those cycles that do not intersect 
the wave's support. The waves move (changing supports in time) in a way that degrees 
are invariant. They are even additive under some circumstances, by next lemma, which 
allows for algebraic manipulations. 

Lemma 12. Let cf>i and <f>2 be two continuous states on X with supports X\ and X2, 
with no two nodes from X\ and X2 being neighbors. Let <j) — (f>i + (f>2 be a state on X, 
then <j) is a continuous state on X, which satisfies h((f>) — h((f>i) + h((f>2). 

Proof. The continuity of <fi inside X\ and X2 is inherited form the continuity of </>i and 
4>2- If x\ £ X\ and X2 ^ X\ are neighbors, then X2 ^ X2- This means <fi(x2) = 0, 
which makes <fi on the pair (x±, X2) is continuous. The same argument works for a pair 
of neighbors in and out of X2- For two neighbors both outside X\ and X2, on which <\> 
is 0, the continuity also holds since the values have to be both 0. Let a = 2~2iLi [ a ii hi] 
be a cycle in X, then aflli and a n X2 are two non neighboring subsets, and: 

K 

h(cl>) = l/nJ2(<t>(h)-<t>(ai)) 

i=l 

For these pairs of neighbors ai and hi, there could be at least one in X\, which sum up 
to be h(<pi), or at least one in X2, which sum up to be h(<p2), otherwise, both are in 
neither X\ or X2, which sum up to 0. Therefore, h{4>) — h(<fii) + h((j>2)- □ 

Corollary 3. If states <f>i, • . • , <f>k have distinct and non-neighboring supports, then 

An important property of the narrow hallways V is has the topological type of a 
planar graph G; specifically, G is a deformation retraction of T>, with retraction map 
r : V -> G and injection map i : G -> V. Suppose Hi(K r (X)) = H x (p) = © g Z, 
and {[ai], . . . , [a g ]} is a basis for Hi(K r (X)), accordingly, {7r*([ai]), . . . , 7r*([a s ])} 
is a basis for Hi(T>). Since the degree of u on a cycle in lZ r (X) is totally determined by 
integers deg(u, ai), . . . , deg(u, a g ) by Lemma|3] we only need to focus on controlling 
the degree on a basis. 



25 



One problem we care about is whether one can realize every possible degree. In 
other words, the question could be reformed as whether the map h is surjective. Specif- 
ically, is it possible to realize a continuous state u, such that deg(u, [a]) = /([a]), 
where / : Hi(TZ r (X)) —> K is any integer valued linear map satisfying f([a + f3]) — 

/([«])+/» 

Our last theorem concerns this ability to program pulses in the network for cus- 
tomizing the response. 

Theorem 8. The map h is surjective: if[f] € H 1 (lZ r (X)), then there exist a continues 
state u on X, such that h(u) = [/]. 

Proof. We start by selecting a specific basis for Hi (G), using the standard basis of the 
complement of a spanning tree T: each remaining edge corresponds with an element 
in a basis of Hi(G). Let this basis be {[a'jj, . . . , [a']}, and the edges in corresponding 
sequence be ei, . . . , e g , where each a is contained in only one element a^. For each i, 
there exist at least one single wave fa that is supported only on a subnetwork in r _1 (e^), 
and satisfies h(fa)([aj]) — §ij, and h(— fa)([(Xi\) = — 1. From the density assumption 
on the network X, those waves can be supported on non-neighboring subnetworks 
and therefore we can sum [/]([«]) of them up to obtain a continuous state </>£ such 
that h((j>' i )([a\) = [/]([o;])^ij by Corallary[3] From the same argument, (/>[ is a 
continuous state which maps to [/] under h. □ 

7 Link Failure Analysis 

Reliability of links is a serious issue for achieving stability of WSN JT] [22) ; in practice, 
stability is not guaranteed, as wireless communication quality is unpredictable under 
different environmental and other physical conditions ll23ll . For our GHM system, it is 
important to keep communication stable, especially the links between sensors of state 
and state 1, since they will determine those nodes' state at the next time step. 

In this section, we will assume that every link works well with a fixed probability 
p s , as a more practical GHM system. By modifying our simulation accordingly, we 
observe that most of the nodes goes to state after the first several steps, as before. Af- 
terwards, either the system dies out if there is no defect (either local or global defects), 
or wavefronts are generated around local defects. But these defects do not guarantee 
the system's periodicity, since link failure might result in their dying, with a probability 
associated with p s . 

For a fixed network X, if given an initial state u with at least one local defect, 
the probability that one local defect dies after T time steps is a function / of X, u, 
T and p s . The smaller p s is, the bigger the probability of defect dying. Meanwhile, 
f(X, u, T, Ps) is an increasing function of T, which approaches 1 as T goes to infinity. 

Although local defects die eventually almost surely, it does not affect pattern prop- 
agation. For a continuous state of waves with no local defect, which are what remain in 
the network after all local defects die, it could either be in a trivial cohomology class, 
which will die out after a while, or has at least one global defect. As in the latter case, 
wave propagation is not necessarily the same as in the deterministic model, since a 
state wavefront may not turn into a state 1 wavefront. However, even this wavefront 



26 



does not update to state 1 as a whole, it is of great chance that at least one of the nodes 
on the wavefront successfully update to state 1 (which still makes a global defect), and 
therefore will gradually correct the neighbors states by contact. 

8 Conclusion 

In this paper, we provide a decentralized, coordinate-free, energy-efficient intruder- 
detection protocol based on the Greenberg-Hastings cyclic cellular automata. The sys- 
tem could easily be adapted to real indoor environment if using sensing devices func- 
tioned with communication and proper sensing ranges. It displays coherence in the 
sense that it is a self-assembling system with random initial conditions; its efficiency 
comes from low power-consuming property inherited from the scheme of the CCA. 
Demonstrations in f|2]are evidence that the system behaves as intended, and this paper 
gives both intuition and rigor about how and why the system works: 

• Wave patterns are explained as a topological phenomenon, determined and de- 
scribed by the existence of defects with nonzero degree. 

• Assigning to wavefronts a cohomology class reveals the qualitative structure of 
the wavefront patterns, greatly clarifying certain classical results about CCA on 
lattices. 

• A non zero restriction of a cohomology class to a subdomain corresponds with a 
set of strategies with which the evader could win the evasion game; meanwhile, 
a zero restriction stands for the failure of the evader: the cohomology class is 
the obstruction for the pursuer to win. This pleasantly resonates with the role of 
cohomology in obstruction theory. 
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